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Abstract 

A  new  conceptualization  '  of  time-frequency  (t-f)  energy  distributions  is 
discussed  in  this  paper.  Many  new  t-f  distributions  with  desirable  properties 
may  now  be  designed  with  relative  ease  by  approaching  the  problem  in  terms  of 
the  ambiguity  plane  representation  of  the  kernel.  Careful  attention  to  the  design 
principles  yields  kernels  which  result  in  high  resolution  t-f  distributions  with  a 
considerable  reduction  of  the  sometimes  troublesome  cross  terms  observed 
when  using  other  distributions  such  as  the  Wigner  Distribution  (WD)  .  When 
these  new  t-f  distributions  are  applied  to  some  common  signals,  fascinating  new 
details  emerge.  Examples  are  provided  for  joint  sounds  and  marine  mammal 
sounds. 

Background 

Signals  of  practical  interest  often  do  not  conform  to  the  requirements  of 
realistic  application  of  Fourier  principles.  The  approach  works  best  when  the 
signal  of  interest  is  composed  of  a  number  of  discrete  frequency  components  so 
that  time  is  not  a  specific  issue,  (e.g.  a  constant  frequency  sinusoid)  or  somewhat 
paradoxically,  when  the  signal  exists  for  a  very  short  time  so  that  its  time  of 
occurrence  is  considered  to  be  known  (e.g.  an  impulse  function).  Much  of  what 
we  are  taught  implies  that  signals  that  can  not  be  satisfactorily  represented  in 
these  ways  are  somehow  suspect  and  must  be  forced  into  the  mold  or 
abandoned. 

It  has  been  quite  difficult  to  handle  nonstationary  signals  such  as  chirps 
satisfactorily  using  conceptualizations  based  on  stationarity.  The  spectrogram 
represents  an  attempt  to  apply  the  Fourier  Transform  for  a  short-time  analysis 
window,  within  which  it  is  hoped  that  the  signal  behaves  reasonably  according 
to  the  requirements  of  stationarity.  By  moving  the  analysis  window  along  the 
signal,  one  hopes  to  track  and  capture  the  variations  of  the  signal  spectrum  as  a 
function  of  time.  The  spectrogram  has  many  useful  properties  including  a  well 
developed  general  theory.  The  spectrogram  often  presents  serious  difficulties 
when  used  to  analyze  rapidly  varying  signals,  however.  If  the  analysis  window 
is  made  short  enough  to  capture  rapid  changes  in  the  signal  it  becomes 
impossible  to  resolve  frequency  components  of  the  signai  -which  are  close  in 
frequency  during  the  analysis  window  duration. 

The  WD  has  many  important  and  interesting  properties.  It  provides  a  high 
resolution  representation  in  time  and  in  frequency  for  a  non  stationary  signal 
such  as  a  chirp.  In  addition,  the  WD  has  the  important  property  of  satisfying  the 
time  and  frequency  marginals.  However,  its  energy  distribution  is  not  non- 
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negative  and  it  often  possesses  severe  cross  terms,  or  interference  terms, 
between  components  in  different  t-f  regions,  potentially  leading  to  confusion  and 
misinterpretation. 

Both  the  spectrogram  and  the  WD  are  members  of  Cohen's  Class  of 
Distributions.  Cohen  has  provided  a  consistent  set  of  definitions  for  a  desirable 
set  of  t-f  distributions  which  has  been  of  great  value  in  guiding  and  clarifying 
efforts  in  this  area  of  research.  Among  the  desirable  properties  are  the  absence 
of  negativity  of  the  energy  values  (the  spectrogram  has  this  property,  the  WD 
does  not)  and  the  property  of  proper  time  and  frequency  marginals  (the  WD  has 
it,  the  spectrogram  does  not).  A  recent  comprehensive  review  by  Cohen  provides 
an  excellent  overview  of  time-frequency  distributions  and  recent  results1. 

Any  member  of  Cohen's  Class  of  distributions  can  be  obtained  by  combining 
a  kernel  with  an  instantaneous  autocorrelation.  In  the  ambiguity  plane.  Some  of 
the  relationships  are  shown  below. 


A(v,  t)g(v,  T) 


S(v,f)  *  y(v,f) 


Usual  computation 
foilowad  by  F 


W(tf)  *  *  t  G(t.O 

W(t,f)=  Jx(t+,c/2)x*(t-T/2)e-j2*fT  dt  (Wigner  Dist.) 


<0* 


A(v,t)=  Jx(t+t/2)x*(t-T/2)e'j25lf'cdx  (Sym.  Ambig.  Fen)  (2) 
*  All  Integrals  are  over  the  entire  real  axis 
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Many  investigators  have  considered  the  WD  to  be  the  basic  t-f  distribution. 
Its  principal  faults,  non-negativity  and  the  cross  terms  have  been  considered  to 
be  facts  of  life  which  can  only  be  dealt  with  by  performing  additional 
computations.  Recent  developments  in  our  group  have  convinced  us  that  there 
are  many  t-f  distributions  with  properties  nearly  as  desirable  as  the  WD,  but 
with  considerably  reduced  interference  terms.  Figure  1  illustrates  the  effect  of 
three  choices  of  kernels,  the  WD  kernel  (Figure  la.),  the  spectrogram  kernel 
(Figure  lb.)  and  a  new  kernel2  designed  to  shape  the  generalized  ambiguity 
function  so  as  to  reduce  cross  terms  (Figure  lc.).  The  kernel  for  this  distribution 
is  of  the  form 

gT  (v ,  t)  =  exp(-v2x2/a)  (3) 

where  o  is  the  parameter  to  control  the  amount  of  interference  suppression.  We 
shall  discuss  distributions  in  this  paper  which  have  kernels  with  similar 
characteristics. 

Reduced  Interference  Distributions 


Figure  Id.  is  the  ambiguity  function  of  two  different  frequency  sinewave 
segments  occuring  at  different  rimes.  The  time  domain  representation  of  these 
two  signals  is  shown  in  front  of  the  t-f  distribution  obtained  assuming  the 
kernels  immediately  above  them  were  multiplied  by  the  ambiguity  function  of 
the  two  signals.  Computation  was  in  the  discrete  form.  Figure  le.  shows  the 
result  obtained  using  the  WD  kernel,  Figure  If.  shows  the  result  obtained  using  a 
possible  spectrogram  kernel  (two-dimensional  gaussian)  and  Figure  lg.  shows 
the  results  using  gr(v,t)  .  One  can  see  that  there  are  strong  interference  terms 
between  the  two  signals  in  the  WD  representation.  The  spectrogram  has  no  cross 
terms,  but  the  t-f  representation  is  not  a  faithful  representation  of  the  time 
duration  and  frequency  spread  of  the  two  signals.  Figure  lg  shows  a  quite 
faithful  representation  of  the  two  signals,  however.  The  rime  and  frequency 
resolution  is  nearly  the  same  as  the  WD,  but  the  interference  terms  are  not 
visible. 

A  descriptive  name  has  been  given  to  distributions  whose  kernels  possess 
the  four  characteristics  which  lead  to  the  much  reduced  interference  terms.  The 
name  Reduced  interference  Djstribution  or  RID  has  been  proposed-.  The  reason 
that  the  RID  approach  is  successful  can  easily  be  seen  by  looking  at  the 
ambiguity  function  in  Figure  Id.  The  desired  components  of  the  ambiguity 
function  fall  near  the  origin.  Th<*  interference  terms  are  far  from  the  origin4. 
The  WD  kernel  includes  them  with  equal  weight  along  with  the  desired  terms 
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because  of  its  unity  value.  The  RID  kernel  suppresses  these  interference  terms, 
however,  and  retains  the  desired  terms. 

Figure  2  represents  a  more  ambitious  comparison  of  the  WD,  spectrogram 
and  RID  kernel  results.  Figures  2a.  2b  and  2c  show  the  results  for  the  WD,  the 
spectrogram  kernel  and  the  RID.  The  signal  is  composed  of  a  chirp  plus  FM 
signal.  The  WD  represents  the  signals  well,  but  strong  interference  terms  are 
observed  between  the  two  signals5.  The  spectrogram  represents  the  chirp  fairly 
well,  but  the  FM  signal  is  very  poorly  represented,  however.  The  RID  result  is 
essentially  as  good  as  the  WD  result  for  the  desired  parts  of  the  t-f  distribution, 
but  with  hardly  any  interference  terms,  in  sharp  contrast  to  the  WD. 

The  required  properties  of  RID,  a  list  of  desirable  properties  and  some 
distribution  comparisons  are  shov,n  in  Tables  1  and  2  the  Essential  Reduced 
Interference  Distribution  Requirements. 


Applications 

Some  comparisons  of  RID  with  spectrograms  and  WDs  are  appropriate. 
Figure  3.  compares  the  Spectrogram  (a),  the  WD  (b)  and  the  RID  (c)  for  2  cycles 
of  a  sinusoid  separated  by  2  cycles,  4  cycles  8  cycles  and  16  cycles  as  one 
progresses  from  left  to  right.  This  is  for  a  256  point  window.  Notice  that  the 
spectrogram  completely  misses  the  on-off  of  the  sinusoidal  pulse.  It  correctly 
reflects  the  harmonic  structure  of  the  on-off  pulsing,  however.  The  WD  reflects 
the  harmonic  structure  in  frequency  and  the  time  structure  as  well,  but  with 
some  breaking  up  due  to  interference  terms.  The  RID  faithfully  reflects  the 
harmonic  structure  and  the  time  structure.  Figure  4  compares  spectrogram, 

WD  and  RID  for  an  off-on  sine  with  increasing  duty  cycle  time  for  window 
lenghts  of  64  128  and  256  points.  The  WD  and  RID  perform  well  here.  TheWD 
interference  terms  are  not  as  evident  due  to  the  finite  analysis  window  lengths. 
The  results  in  Figures  3  and  4  were  motivated  by  Watkin's  paper6  on  limitations 
of  the  spectrogram  and  include  recreations  of  his  spectrogram  results. 

Figure  5.  compares  spectrogram  and  RID  results  for  a  TMJ  click  during  jaw 
movement.  Figure  6.  shows  RID  results  for  a  bottlenose  dolphin.  Note  the  relative 
invariance  of  the  t-f  patterns  compared  to  the  time  series  underneath.  The 
spectrum  gives  little  hint  of  the  t-f  structure.  Figure  7  shows  RID  results  for 
spermwhale  clicks.  Note  the  repeated  triplets  in  the  data  from  the  two 
hydrophones. 

Conclusions 
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1.  RID  produces  very  useful  and  interpretable  results  for  natural  signals.  Unique 
signatures  for  underwater  sounds  seem  to  be  an  attractive  feature  and 
capability  of  RID. 

2.  RID  well  represents  what  one  would  conceptually  expect  from  an  ideal  time- 
frequency  distribution  in  many  situations. 

3.  The  precise  time  and  frequency  structure  of  RID  may  allow  source  and 
multipath  structures  to  be  disassociated. 

4.  Care  must  be  taken  in  proper  computation,  display  and  interpretation  of 
results. 

5.  There  are  conditions  where  the  performance  of  RID  is  not  as  good  as  a 
specialized  distribution  or  approach  designed  to  the  signal  and  situation. 

6.  The  RID  possesses  almost  all  of  the  desirable  properties  of  a  time- 
frequency  distribution  and  is  particularly  useful  in  multicomponent  situations 
where  the  reduced  interference  property  is  very  desirable. 

7.  The  Wigner  distribution  still  has  an  important  place  in  theory  and  in 
application,  particularly  when  interference  is  not  an  issue,  when  one  desires  to 
use  time-varying  filtering  and  synthesis  techniques  and  when  signal  detection  is 
an  issue.  It  may  also  be  desired  if  cross  terms  are  useful. 
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'PO.  ncnnegativity  :  p-(t, /;  g)  >  0  vt, } 

QO.  g(v,r)  is  the  ambiguity  function  of  some  function  vu(t). 

PI.  realness  :  pz(t,f:g)€R 
Ql-  $(«.')  =  S*C-v.  -r) 

P‘2,  time  shift  :  s(t)  =  z(t  -  t0)  =>  ps(t./:j)  =  />:(£  -  tQ,f:g ) 

Q2.  g(v,r)  does  not  depend  on  t. 

P3.  frequency  shift  :  s(0  =  =>  /;  g)  =  p2(t.  /  -  /o : 

Q3.  g(v,r)  does  not  depend  on  /. 

P4.  time  marginal  :  J"  Wz(t,f)df  =  r(t)**(0 
Q4.  5(u,0)  =  1  Vv 

PS.  frequency  marginal  :  j* Pz(t,f:g)dt  =  Z(f)Z*(f ) 

Q5.  5(0, r)  =  1  Vr 

f  fpz(t,f:g)df 

_ ...  ,  _  .  J  _  .r  /^S 


P6.  instantaneous  frequency  : 


Q6.  Q4  and 


3s(v,r) 


J  Pz(t,f;g)df 


=  /,(£) 


r=0  =  ov” 


/  tpz(tj:g)dt 

P 7.  group  delay  :  - -  =  t$(/) 

Q7 .  Q5  and  - - - |u_0  =  0  Vr 

dv 

P8.  time  support  :  z(t)  =  0  for  |t|  >  te  =>  Pz(t,f'.g)  =  0  for  |t|  >  tc 
Q8.  ^(t,r)  =  j g(v,T)e-i2levtdv  =  0  for  |r|  <  2|t| 

P9.  frequency  support  :  2(/)  =  0  for  |/|  >  fc  =>  pz(t,f:g )  =  0  for  |/|  >  /c 
Q9.  J g(v,r)e^TfTdr  =  0  for  |v|  <  2|/| 

P10.  Reduced  Interference 

Q10.  j(v,r)  is  a  2-0  low  pass  filter  type. 


Tabic  I:  Distribution  properties  and  associated  Kernel  requirements. 


Distribution 


PO  PI  P 2  P3  P4  PS  P6  P7  P8  P9  P10 


Wigner 

Rihaczek 


Re{Rihaczek}  cos(irvr)  x  x  x  x 

Exponential  (ED)  e-v  T/a  xxxxxxx 

Spectrogram  Xw(v,t)  of  a  window  «(t)  x  x  x  x  x 

sin(#T/r) 

p-> _ l  _ _ J _ _  u  u  w  V 


xxxxxxxxx 


X  X  X  X  XX 

xxxxxxxxx 


Born- Jordan* 


xxxxxxxxxx 


•j 

Windowed-ED“  «“v !c  *  W(v)|v=u 


XXXXXXXXXX 


*  belongs  to  the  RID. 


Table  2:  Comparison  among  several  distributions. 


RID  kernel  g(v,T) 


RID  kernel  Requirements  for  Time  and  Frequency  Support 

Essential  Reduced  Interference  Distribution  Requirements 


FIGURE  3.  A  simulated  sinusoid  is  pulsed  2  cycles  on,  2  cycles  off;  2  cycles 
on,  4  cycles  off;  2  cycles  on,  8  cycles  off;  2  cycles  on,  16  cycles  off  for  the 
(a)  spectrogram,  (b)  Wigner  distribution,  and  (c)  RID.  The  analysis  window 
length  is  236  points  in  all  cases. 


RID:  WINDOW  l  ENOTIUM  POINTS 


RIO:  WINDOW  LENGTH-128  POINTS 


RID:  WINDOW  LENGTIUZ56  POINTS 


time  -> 


